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Magnetohydrodynamic Turbulent Cascade of Coronal Loop Magnetic Fields 
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The Parker model for coronal heating is investigated through a high resolution simulation. An 
inertial range is resolved where fluctuating magnetic energy Ej^j{k±) oc fcj^'^ exceeds kinetic energy 
E^(k_i_) oc fc^ ■ . Increments scale as 5h^ ~ l~ ' and 5u^ ~ t^ ' with velocity increasing at small 
scales, indicating that magnetic reconnection plays a prime role in this turbulent system. We show 
that spectral energy transport is akin to standard magnetohydrodynamic (MHD) turbulence even 
for a system of reconnecting current sheets sustained by the boundary. In this new MHD turbulent 
cascade, kinetic energy flows are negligible while cross-field flows are enhanced, and through a series 
of "reflections" between the two fields, cascade more than half of the total spectral energy flow. 

PACS numbers: 47.27.Ak, 47.27.ck, 96.60.pf, 96.60Q- 



The heating of solar and stellar atmospheres is an out- 
standing astrophysical problem [1| . The solar corona has 
temperatures (> 10^ K) up to three orders of magni- 
tude higher than the underlying photospheric and chro- 
niospheric layers, sustained by an energy flux of about 
10^ erg cm~^ s~^ Q. 

Most of the coronal x-ray and extreme ultraviolet 
(EUV) radiation is emitted in loops, bright structures 
threaded by a strong axial magnetic field connecting 
photospheric regions of opposite polarity. In the sce- 
nario proposed by Parker [3|, magnetic field lines are 
braided by convective photospheric motions that shuffle 
their footpoints, leading to the "spontaneous" develop- 
ment of small-scale current sheets where the plasma is 
heated. 

It has long been proposed Q that this scenario can 
be regarded as a magnetohydrodynamic (MHD) turbu- 
lence problem as photospheric motions stir the magnetic 
field lines' footpoints, and these motions are transmit- 
ted inside by the field-line tension stirring in this way 
(anisotropically) the whole plasma akin to a body force. 
Simulations have indeed revealed that the system ex- 
hibits many properties of an authentic MHD turbulent 
system, including the formation of field-aligned current 
sheets, power-law spectra for the energies, and power-law 
distributions for energy releasejjeak dissipation and du- 
ration of dissipative events [5|, |6| . Furthermore in recent 
papers [6| we have developed a phenomenological scaling 
model for this turbulent cascade where at the large scales 
nonlinearity is weak (i.e., depleted similarly to |7|) and 
at the small scales strong Q . 

However this new, line-tied, turbulent regime is quite 
distinct from the classical MHD turbulence system where 
energies are in approximate equipartition, as here fiuc- 
tuating magnetic energy dominates over kinetic energy 
throughout the inertial range. 

Certainly, this system can also be seen as a set of re- 
connecting current sheets sustained by the boundaries, 
and a large fraction of the kinetic energy might be con- 
tributed by the magnetic field through reconnection itself 
rather than from cascading large-scale kinetic energy. 



It is therefore crucial to understand if turbulence is an 
appropriate framework to model this problem, namely, 
can a set of reconnecting current sheets be described 
in terms of turbulence? The influence of turbulence on 
magnetic reconnection is an active research topic (e.g., 
see [9| for a model), but are energy fluxes in a system 
where magnetic reconnection plays a prime role similar 
to those of MHD turbulence? We present here a novel 
investigation of the energy flows between different scales 
and fields in Parker's model in order to determine how 
different the spectral fluxes are and what properties they 
share with the standard MHD turbulence case Il0l-[l2i . 



A coronal loop is modeled in Cartesian geometry as a 
plasma with uniform density po embedded in a strong ax- 
ial magnetic field Bq directed along z (see p for a more 
detailed description of the model and numerical code). 
Magnetic field lines are line tied at the top and bottom 
plates where a large-scale velocity field is imposed. In 
the perpendicular direction {x-y) periodic boundary con- 
ditions are used. The dynamics are modeled with the 
(non-dimensional) equations of reduced magnetohydro- 
dynamics (RMHD) [ll]: 



dtn + (u • VJ u = -V^P + (b • VJ b + c^d.h + -V>, 

H 

dth + (u • vj b = (b • vj u + c^a.u + 4v>, (1) 

ri 



with V^- u = V^- b = 0. Here gradient and Laplacian 
operators have only orthogonal {x-y) components as do 
velocity and magnetic field vectors {uz=hz=Q), while P is 
the total (plasma plus magnetic) pressure, c^ is the ratio 
between the Alfvcn velocity of the axial field {Bq / ^J At: p^) 
and the rms of photospheric velocity (1 kms~^) and R is 
the Reynolds number. In the simulation presented here, 
c^ = 400, R = 2000, and the domain spans < x, y < £, 
< z < L, with £ = 1,L = 10 and 1024^ x 512 
grid points. Given the orthogonal Fourier transform 
u(x^,z) = Yj\^ u(k^ , z) e^'^i'^i with k^ = 27rn^/^ and 
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FIG. 1. Magnetic (-Ej^j) and kinetic (-B^) energy spectra as 
a function of the orthogonal wave number n^ (k^ — 2im_^^ ) . 
Inset: Increments, velocity increasing at high wave numbers. 
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FIG. 2. (Color onhne) Spectral fluxes [Eqs. ((6|, Q, (0, and 
(|10|) ] showing incoming and outgoing energy transfers between 
different shells and fields to/from shell Q=20. 



i\ G 1?, the shell- filtered field u^^. is defined as [Il| 

u,(x,,z)= Y. fl(k,,^)e^'^-''-, (2) 

i.e., it has only components in the "shell" K with wave 
numbers K-l< jrij^l <K. The boundary photospheric ve- 
locities at z = 0,L are given random amplitudes for all 
wave numbers 3<n^ <4 and then normalized so that the 
rms value is 1/2 [seeQ. As a result the forcing boundary 
velocity has only components in shells 3 and 4. 

Since fields filtered in different shells are orthogonal, 
and indicating the volume integrals with (...), the equa- 
tions for kinetic and magnetic energies £'(^' — (u^ /2), 
E[K) 



;b2 /2) in shell K follow from Eqs. (P): 
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'' b / J ub bb 






(3) 



(4) 



These are obtained as for the three-periodic case [11| . 
except for terms ex dz in Eqs. ([T]) that contribute the 
integrated Poynting flux S"'^^ (the work done by photo- 
spheric motions on magnetic field lines' footpoints) en- 
tering the system at the boundaries in shell K: 



g(K)^ 



da u,^- b — 



z=L 



da u ^ • b 



z=0 



(5) 



This docs not cancel out along the nonperiodic axial di- 
rection z. As photospheric velocities have only compo- 
nents in shells 3 and 4 S''^-' vanishes outside these two 
shells {the injection scale). In similar fashion the dissi- 
pative terms V[^^ = (|Vu,J^), V['^'> = (|Vb^,|^) con- 
tribute only at dissipative scales with large K. 

Between the injection and dissipative scales only the 



following terms contribute; 

ri^'^) = -(u,,.(u.v)uj 

1 



Ai^^^c,{u,-dzh,)-^S^'<^ 
7^(Q,i^)^_(b^,.(u.v)bJ 



7-(Q 



K) 



:b,-(b-V)u, 



Vl(^'5 



Q,K 



A[P=cJh,-dzn,. 



2^ 



(6) 
(7) 
(8) 
(9) 

(10) 
(11) 



These represent energy fiuxes between fields in different 
shells. In fact given two fields v and w (either velocity 
or magnetic fields) the relation T"'-'^'^^-' = —'J'(^'Q) holds 
for the flux between shells Q and K, and together with 
Eqs. (©-(Ill) define the fluxes [III- For example, T}^'^^ 
represents conversion of kinetic energy in shell Q to mag- 
netic energy in shell K. 

A^^'' is the flux due to the linear terms ex d, in 
Eqs. ([IJ. It does not transfer energy between differ- 
ent shells, but only in the same shell between fields u 
and b. Notice that it has been symmetrized so that 
A[^^ = —Alff\ as can be verified integrating by parts. 

Equations ([31)- jl]) show that at the injection scale 
(shells 3 and 4) the photospheric forcing supplies en- 
ergy at the same rate to both kinetic and magnetic ener- 
gies, i.e., the forcing injects Alfven waves unlike standard 
forced MHD turbulence where a mechanical force injects 
only kinetic energy. 

The simulation is started with vanishing orthogonal 
velocity and magnetic fields (u = b = 0) and a uni- 
form axial magnetic field B = i?o Gz inside the compu- 
tational box. As shown in our previous works [6| the 
constant forcing velocity at the boundary advects mag- 
netic field lines generating a perpendicular component 
b. This initially grows linearly in time before saturating 
nonlincaritics develop, and kinetic and magnetic energies 
then fluctuate around a mean value, with magnetic field 
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FIG. 3. (Color online) Energy flows between scales (log 
bands) of the magnetic field. Log bands (pft] = (2'''~^, 2*"^]; 
see text. 
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FIG. 4. (Color online) Energy fluxes from magnetic to ve- 
locity fields scales. Log bands, e.g., (pa] = [3, 4] and (p6] = [17, 
32] 



fluctuations dominating; E^,^ /E^ ^ 87 for the simulation 
presented here. 

The magnetic and kinetic energy imbalance is reflected 
in the energy spectra shown in Fig. [T] Both spectra are 
peaked at the injection wave numbers 3 and 4, but be- 
yond n^= 5 an inertial range is resolved where both spec- 
tra exhibit a power-law behavior with a steep index for 
the magnetic energy (—2.7) and a flatter one for the ki- 
netic energy (—0.6). Since increments are obtained from 
band- integrated spectra [e.g., 6b'^{kJ ~ k^EM{kJ], this 
implies 



Sk 



0.85, CT - -0.2, (12) 



i.e., while magnetic energy decreases at small scales ki- 
netic energy increases as it is expected to do for magnetic 
reconnection [9[ where velocity is concentrated at small- 
scale current sheets [a]. 

Energy fluxes are shown in Fig. [2] As their behavior 
is similar along the whole inertial range, here we plot the 
fluxes in and out of shell Q = 20. With respect to the 
equipartition case (EQPT) [11| , the most striking differ- 
ence is the small value of the transfers between kinetic 
energy shells 7~(Q=20,/-f)^ negligible with respect to the 
others. Indeed, the velocity eddies are not distorted by 
other velocity eddies as they are too weak compared to 
the strength of the magnetic field: both the orthogonal 
component b and the strong axial field Bo are responsible 
for shaping the velocity field. 

We analyze first the energy fiows between shells of the 
magnetic field T^W^^o.if)^ They are negative for all K<Q 
and positive for all K>Q meaning that the field is receiv- 
ing energy from shells at smaller K and giving energy to 
shells of greater K. In contrast to EQPT for K<Q there 
is an almost constant small contribution from smaller 
K shells. As in [ll[ a similar "tail" is also present for 
T^^^'^'^'^\ the magnetic field in shell 20 is receiving en- 
ergy from smaller K shells of the velocity field and trans- 
ferring it to larger K shells. T}^=^"'^^ = -T}f'^=^"^ 
has the corresponding behavior. 



The large peaks at K=20 represent the conversion of 
magnetic to kinetic energy in the same shell and are due 
to A'-^=^°^ > {A[^^=-A[^^). Its large value is linked 
to the field-line tension of the dominant axial field Bo as 
A^^ is obtained from the linear terms oc c^dz in eqs. p]). 
this is the Alfvcn propagation term that in presence of 
a b^ contributes with a velocity u^ of the same shape. 
In Fig. [2] the values of the cross-field fluxes 7^„ and T^^ 
for K=20 without these contributions are shown with a 
diamond and a solid circle, respectively 

The "tails" shown in Fig. [2] are also present for higher 
values of Q. In [11|, they were present only in cross-fields 
transfers, while here a comparable tail appears also in 7^j_ . 
The tail in this work is due to the steeper spectrum of 
the magnetic field; this induces a higher spectral trans- 
fer at low K because of the higher value of b^ in flux 
X, [Eq. (HI)] with respect to the EQPT case ^. This 
feature has been indicated as evidence of the nonlocal na- 
ture of energy transfers in MHD turbulence [11| because 
the cumulative transfers of farther shells seem to be more 
important than those of closer shells. 

However, in scaling models of MHD turbulence [3, Q 
and of Parker's model [g] scales are defined as log bands 
of shells. In fact a single shell does not represent a 
scale since, from the uncertainty principle, its associated 
field [Eq. ([2])] is delocalized in space and cannot represent 
an eddy, the building block of K41 phenomenology. 

Thus in order to understand how energy fiows across 
scales we must use such log bands [l2| . Log band {p] is de- 
fined as the shells included in (p/2,p], equally spaced on a 
logarithmic scale. Considering p„ = 2"~^ we will indicate 
these intervals with their index n: (p„] = (2"~^,2"~^] 
(n=l,...,10). With a 1024^ grid we have "only" 10 dis- 
tinct intervals. Figures [3] and S] show the fluxes summed 
over these log bands of shells, e.g., T^^''^^ =j-ap,]XPk]) 
for q=6, 7 and 8. The injection scale is now (p3] = [3,4], 
while (p6] = [17,32], (p7]-[33,64], (p8] = [65,128], etc.. 

Figures [5]|3] show that the apparently dominating con- 
tributions of distant shells (tails in Fig. [J) strongly de- 



E. 



E„ 



k 



(log bands) 



-? ? r ? ? ■' k 



..Y ••■>.■■■ ••■>.■■■ •••>.■■■ ••■>.■ 

FIG. 5. Diagram of energy flows between fields and scales. 



crease when considering log bands. In these bands 
the number of shells increases exponentially at higher 
wave numbers and the aggregate effect of local transfers 
asymptotically dominates J12l . [ij] . Recently an analyti- 
cal upper bound has been set for the locality of energy 
transfers [12|, although within these bounds the energy 
transfers can be quite spread out 15| . In fact while cross 
energy transfers (Fig. |4]) arc quite local as energy flows 
between neighboring scales decrease swiftly, the transfers 
between scales of the magnetic field (Fig. [3|) are more 
spread out. However for the Parker problem we do not 
observe a direct flow of energy between the forcing scale 
and the small scales UM. 



Overall energy is transferred from larger to smaller 
scales in a similar fashion as in MHD turbulence with en- 
ergies in equipartition, except for the velocity field that is 
too weak compared to the magnetic field. As a result the 
magnetic field creates and shapes the velocity field. In 
fact from cross-field flows (Fig. 2]) we sec that magnetic 
field line tension enhanced by line tying, and predomi- 
nantly represented by the fluxes A^^^^ , converts magnetic 
energy to kinetic energy at the same scale. In turn kinetic 
energy at larger scales is converted to magnetic energy at 
smaller scales, due to the magnetic stretching term. The 
magnetic advection term transfers a (smaller) fraction of 
energy toward smaller magnetic field scales (Fig. [3]). A 
pictorial summary of the cascade is shown in Fig. [5] (mag- 
netic flux spread not shown); the repeated conversion of 
kinetic to magnetic energies by the cross-field fiows effec- 
tively cascades energy toward the small scales. 

The upper bound for the locality of energy transfers 
found by [l3] has been restricted to the case < a; ^ < 1 
in analogy to the hydrodynamic case [iJ] but this condi- 
tion is overrestrictive for MHD and we show that those 
bounds are valid also for the case [Eq. ([12])] where a^ ^ 
0.85, but a^ ^ —0.2 instead of the standard a^ ^^ ~ 1/3. 

If [K] is the log band of shells [K/2, K], heuristic scal- 
ings [Eq. (|12p ] can be written more precisely for the 
generic band-summed field a^^'l = ^ a^^ with K € [K] 



[Eq. (ED] as 
/|aM|3)i/3^/^-. 



[Val^'I [3)1/3 ^^-i--„. (13) 



Since the scaling for the derivative is valid independently 
of the sign of o; , following [l^ the bound for the locality 
of cross-field transfers between band-summed fields is 



bxL 

r([Q]:M) 



< ^-(a„+a6)gi-a,^ for Q<A'(14) 



< K 



l-O-u 



'2(76 



for g>A' (15) 



At fixed K, contributions from smaller bands [Eq. ()14p ] 
are negligible if ct^ < 1 so as from bigger bands [Eq. (|15p ] 
for Oj > 0. 

In similar fashion, we obtain the following bound for 
magnetic fluxes: 



r([Q],[^]; 



y([Q],[^]) 



^ ^-('yu+'yt) Qi-'yt ^ for Q < A' (16) 



<K' 



i-o"b Q-(<yu+<yb) 



for Q:$>K (17) 



The requirement for asymptotic locality is still o; < 1 for 
Eq. ([H]), but CT, + cTj, > for Eq. dTT]), ah satisfied in our 
case. 

The sign of the exponents containing a^ remains un- 
altered with respect to the classic case a; „ ~ 1/3, but 
while in Eqs. (|14p and (fT6|) /v ^0-6^ decreases for large 
K in Eqs. (fT5|) and (fT7|) we have a positive power {K^'^ 
and A'^i^, respectively) as in the standard case [l2| with 
K^l^ . As in the hydrodynamic case [IJj the origin of this 
pathological behavior stems from the use of Holder in- 
equality for fiuxes [Eqs. ([6])-([9])] to set the upper bound, 
since only the absolute values of their terms arc con- 
sidered, and any cancellation effect due to the scalar 
products in Eqs. (I6])-(I9|) is neglected. This conclusion 
is reinforced by the fact that for Parker problem fluxes, 
Eqs. (©-([HI) exhibit scalings (not shown) in Q well below 
upper bounds [Eqs. (fT4)) - (fT7)) ]. 
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